Collective oscillations of a trapped Fermi gas near the unitary limit 
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We calculate the oscillation frequencies of trapped Fermi condensate with particular emphasis 
on the equation of state of the interacting Fermi system. We confirm Stringari's finding that the 
frequencies are independent of the interaction in the unitary limit, and we extend the theory away 
from that limit, where the interaction does affect the frequencies of the compressional modes only. 
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The remarkable advances in producing and measuring 
properties of atomic condensates give a strong impetus 
to develop theory to meet the challenges of interpreting 
the experiments. In the case of fermion condensates, one 
is now at the early stages of coming to an understanding 
of the first experimental results [ij, y| • One of the char- 
acteristic properties is the frequencies of normal modes 
of vibration. Stringari Q has developed the theory at 
the unitary limit, finding that the oscillation frequencies 
are independent of the details of the interaction. Here 
we extend the theory away from the unitary limit where 
the interaction has some effect. 

The "unitary limit" is a term to describe a two- 
component Fermi gas with a short-range interaction, 
characterized by a scattering length that is large com- 
pared to the length scale set by the particle density. This 
limit was discussed in 1999, when one of us (GFB) for- 
mulated as a challenge to many-body theorists to clarify 
the structure of the ground state of a fictitious neutron 
matter, interacting with an infinite scattering length 
At the time when the challenge has been issued it was not 
really clear even if such rnatter is stable in principle, as 
either a system of_bosons or a system of three or more 
fermion species 6] is known to be unstable, due to the 
well known Efimov effect. From simple dimensional ar- 
guments it was clear that the energy per particle of a two 
fermion species should be proportional to that of the free 
Fermi gas, with a universal constant of proportionality. 
This may be expressed as 
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where kp is the Fermi wave vector and ^ some pure num- 
ber. Even though estimates of this number have been 
generated at the time 0, the most convincing argument 
that this number is positive indeed came only recently, 
both in theory namely ^ w 0.44, and experiments 
0. Experiments indicate so far that ^ w 0.5 with very 
large error bars. The theory also provides strong argu- 
ments that such a fermionic system should be superfluid 
at zero temperature, while the experiment is still lacking 
in this respect. The information inferred from various 
experiments so far is, at best, of a rather indirect na- 



ture and no direct unambiguous evidence of superfiuidity 
has been reported until now. Very likely, only the direct 
observation of vortices will provide the ultimate experi- 
mental proof that such systems can sustain a superflow. 
There are strong theoretical arguments that such vor- 
tices should be almost as easily visible in the Fermi dilute 
atomic clouds JXLj as in Bose dilute atomic clouds 

In the Green Function Monte Carlo calculations of the 
Refs. 0, ^1 it was established that the ground state en- 
ergy per particle of the homogeneous phase of a system 
of two fermion species, interacting with a very large scat- 
tering length (much larger than the radius of the interac- 
tion ro (namely \a\ 3> tq), in the dilute regime (nrg ^ 1, 
where n ^ + ni — kp/Sn'^) is 
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where F{x) is a universal function and the constants ^ ~ 
0.44, C « 1 and v ^ 1 (the last two values have been 
extracted by us from the numerical results provided by 
the authors We shall discuss the normal modes of 

oscillation for a Fermi system in a harmonic trap, whose 
energy per particle is given by Eq. We note that the 
last term is independent of density and can be dropped 
from the present analysis. 

We assume throughout that the system behaves hydro- 
dynamically, i.e. that the pressure tensor is isotropic. If 
the system is superfluid, then as long as the oscillation 
frequency is below the gap frequency needed to break-up 
a Cooper pair, this condition is expected to be fulfilled. 
The ground state density satisfies the equation 
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where uq is the ground state density, U is the trapping 
potential, and P is the pressure. P is related to the 
energy per particle e{n) by 
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The slow and small-amplitude normal modes satisfy the 
following equation 
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where rii is the oscillating density and 
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gives the speed of sound c in a uniform condensate at 
density n. It is st raig htforward to show that Kohn's gen- 
eralized theorem (stating that the frequencies of the 
dipole modes are the trap frequencies) is satisfied. It is 
convenient to make a change of variable in this equation, 
replacing ni by the variable /i defined by 



"1 = nofi/c^ino). 
After some algebra making use of Eq. (j2Jl we find 
V-(noV/i) = 
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This equation has the formal advantage of being Hermi- 
tian, and thus easier to study in perturbation theory. 

As noted before 0|, the equations of motion admit 
simple scaling solutions for gases in harmonic traps that 
satisfy polytropic equations of state. Before discussing 
a perturbative treatment, we shall present results of an 
analysis under the assumption of a polytropic equation 
of state. As we were completing this work, we learned of 
a similar analysis by Heiselberg, who presents a general 
solution in terms of hypergeometric functions (which re- 
duce to polynomials in this case) The polytropic 
equation of state is 
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where the constant 7 is the adiabatic index of the gas. 
The ground state density is given by 
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Here we introduced 
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where /ip is the chemical potential and the harmonic trap 
potential is given by 
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The corresponding local sound speed is given by 
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TABLE I: Results for a polytropic gas. For A 1 only leading 
terms are shown and £1,2,3 are some constants. 



trap type mode 


/i 




sphericcil L = 1 
A = 1 L = 2 
L = 


X, y, z 
xy, etc. 
1 - 2r^ 


1 
1 

2 

37-1 


axial M = ±2 
A < 1 A/ = ±1 
radial 
axial 


xy, x^ - 
xz, yz 
x^ + ciX^z'^ +C2 

Z^+C3 


2 

1 + A^ 
27 

A2(37-1)/7 



Using these expressions, Eq. 
form 



^ can be written in the 
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Due to the particular polynomial form of fi{r) the eigen- 
functions /i have a polynomial character as well. From 
Eq. ||2Jl one can extract an effective adiabatic index for 
a Fermi gas in the vicinity of a Feshbach resonance 
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Note that the quadrupole mode in the spherical con- 
densate and the transverse quadrupole modes in the de- 
formed condensate do not depend on the equation of 
state. That is because the flow in these modes is incom- 
pressible and the internal energy does not change during 
the oscillation cycle. The frequencies of the monopole 
and of the two compressional modes in the deformed con- 
densate have a dependence on 7, and we can use that to 
estimate the frequency shift. 

Eq. (|15|l shows that the effective adiabatic index is 
larger than 5/3 on the BEC side of the Feshbach reso- 
nance (when a > 0). This behavior implies that the fre- 
quency of the radial oscillations should increase as well 
when going from the BCS to the BEC side of the Fes- 
hbach resonance. This conclusion agrees with the con- 
jecture made by Stringari except that now the fre- 
quency shift has been evaluated explicitly in terms of 
the properties of the system. In the unitary limit, when 
7 = 5/3, these results also agree with previous results for 
the non-interacting Fermi systems in traps |l6| and the 
results for a superfluid Fermi system away from a Fesh- 
bach resonance in a spherical trap The quadrupole 
frequencies obtained using scaling solutions jlj| and the 
sum- rule approach |l8l |. in the limit of a non- interacting 
Fermi gas, are different, and indeed correspond to the 
diabatic limit or coUisionless regime 0|. In this limit 
the sphericity of the Fermi surface is lost during oscilla- 
tions and the cloud behaves like a normal Fermi gas in 
Landau's zero sound regime. 

The polytropic analysis is useful to show the basic de- 
pendence on the system parameters, but the parameter 
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kp is ill-determined, due to the non-uniformity of the 
condensate. We shall therefore use perturbation theory 
to make a more quantitative analysis. Since Eq. JSJ is 
Hermitian, the following formula is variational in the fi , 
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This implies that one can ignore the perturbation in /i 
in calculating the first-order shift in the frequency. That 
shift is then given by 

5u.\_ J\Vh\Hn<fr f f?6 {^) d^r 



To apply this formula, we take the unperturbed densities 
no and ng/c^ from Eq. (10) and (13), setting 7 — 5/3. 
The needed functions are 
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We now have to determine the change in n and n/c^ 
induced by the second term in the energy function Eq. 

Including that term, the chemical potential satisfies 
the equation 
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We may use this expression to evaluate the first-order 
change in ground state number density no- Since the 
frequency does not depend on ^,o^ we may hold it fixed 
in doing the variation. The result is 
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In the same way we include the perturbation in the for- 
mula for and obtain that the corresponding first order 
change in n/c^ is given by 
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We next evaluate the integrals in Eq. H17() . It is 
convenient to change lengths to a dimensionless form, 
scaling them by the transverse condensate radius R = 

a/2/Xo/toWo; 



R{x,y,z) = {x,y,\z). 
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In terms of the scaled variable jl = 1 — f ^ , Eq. I|16(l for 
the unperturbed frequency may be expressed 
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Eq. (12) for the frequency shift becomes 



5uj^ _ C h 
up' ^1/2 miooRa, 



1 



where the dimensionless factor K is given by 
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and kp{0) is the value of the local Fermi momentum at 
the center of the trap. The prefactor in Eq. H24|) dis- 
plays the scaling of the frequency shift with respect to the 
physical parameters of the condensate. As expected, the 
shift is inversely proportional to the combination kpa. 
Finally, it has a nontrivial dependence on ^ and Ci the 
universal parameters defining the energy per particle in 
the vicinity of a Feshbach resonance. 

All the needed radial integrals have the form 
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after integration over angular variables. Since in Eq. I|25() 
both denominators and numerators have the same angu- 
lar dependence, the specific values of the angular inte- 
grals cancels out in the case of incompressible modes. 
We present the results for the various cases in Table 
II. One sees that the shift of the dipole mode vanishes, 
as required by the generalized Kohn's theorem. The 
shift also vanishes for the pure quadrupole modes, for 
reasons noted earlier. The cases of most interest are 
the monopole mode in spherical traps (A = 1) and the 
M = modes in axially deformed (essentially cylindri- 
cal) traps with A ^ 1. The K factors are non- vanishing 
in these cases, but they are rather small, for example, 
Kradiai ~ 0.12. This has the same order of magnitude as 
the factor 1/10 in Eq. (15) for the change in the effective 
adiabatic index. We also note that the sum of factors 
K determining the shifts for the radial and axial modes 
equals the factor K for the pure monopole mode in the 
spherical case. 

The two experimental results available so far 0, |3| are 
still in noticeable disagreement with each other to permit 
a detailed comparison with theory. Nevertheless, both 
experiments show distinctly a qualitative agreement with 
theory as far as the character of the frequency shift is con- 
cerned, in the vicinity of the Feshbach resonance. The 
fact that both experiments seem to favor the adiabatic 
character of the oscillations should not be interpreted 
yet as a confirmation of the existence of superfluidity in 
these systems, since the sphericity of the Fermi surface 
can be maintained by collisions. In this respect the sit- 
uation here is to some extent similar to the expansion 
of a cold Fermi gas |^ . It is also important to de- 
termine experimentally the frequencies of the transversal 
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TABLE II: Results for K. 



trfip type 


mode 




UJ 




spherical 


dipole 


2 


LOO 





A = i 


monopole 


1 - 2r^ 




256 
5257r 




quadrupole 









axial 


M = ±2 


xy, - 







A< 1 


M = ±l 


xz, yz 


LOO 







radial 


x'+y^ + iX^z^- 




1024 




26257r 




axial 


1 - eA^z^ 


Aco 


256 




26257r 



quadrupole modes, since a shift in these frequencies can 
point to a complex structure of the cloud, similar to that 
discussed in Ref. |23 |. 
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Note added. After submitting this manuscript we 
learned of a few another studies of Fermi systems using 
the polytropic equation of state ^] . 
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